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Abstract —How would a cellular network designed for high 
energy efficiency look like? To answer this fundamental ques¬ 
tion, we model cellular networks using stochastic geometry and 
optimize the energy efficiency with respect to the density of 
hase stations, the number of antennas and users per cell, the 
transmit power levels, and the pilot reuse. The highest efficiency 
is neither achieved hy a pure small-ceU approach, nor hy a pure 
massive MIMO solution. Interestingly, it is the combination of 
these approaches that provides the highest energy efficiency; small 
cells contributes by reducing the propagation losses while massive 
MIMO enables multiplexing of users with controlled interference. 

I. Introduction 

Two key goals for the fifth generation (5G) cellular net¬ 
works are improved spectral efficiency (SE) and higher energy 
efficiency (EE) Q. These performance metrics are coupled and 
cannot be treated separately in the design of future networks 
Q. The key to improve the SE and EE is higher spatial reuse; 
that is, more parallel transmissions per km^. There are two 
main densification approaches: 1) smaller cell radius Q and 
2) massive MIMO (multiple input, multiple output) technology 
0. The purpose of this paper is to show that these approaches 
are fundamentally non-competing; in fact, both are needed to 
make future wireless networks truly energy efficient. 

The EE is defined as the benefit-cost ratio of the network: 

Area Spectral Efficiency [bit/symbol/km^] 

!E!E — 2 -1 

Transmit + Circuit Power per Area [Joule/symbol/km J 

This definition reveals that there are three main factors that can 
be modified when optimizing the EE of a cellular network. 

Smaller cells are obtained by densifying the base station 
(BS) deployment. This has a positive effect on the EE in terms 
of increasing the SE and reducing the transmit power (due to 
lower propagation losses). However, the negative effect is the 
increased circuit power due the larger amount of hardware 
infrastructure. Small-cell deployments need to be modeled 
in an asymmetric fashion since the users are non-unifotmly 
distributed over the coverage area. This makes Poisson point 
processes (PPPs) suitable and tractable analytic models 0. 

In contrast, massive MIMO evolves the conventional BS 
technology by replacing the bulky high-gain vertical antennas 
with arrays comprising hundreds of small dipole antennas. 
Contrary to what the term “massive” suggests, massive MIMO 
arrays are rather compact; 160 dual-polarized antennas at 3.7 
GHz fit into the form factor of a flat-screen television By 
processing the antennas coherently, the array can receive and 
transmit a multitude of signals for users at different spatial 
locations. The positive effects on the EE comes from higher 
SE and less transmit power (due to an array gain), but the 
circuit power increases with the number of BS antennas. 
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Fig. I. Illustration of one realization of BS positions from the PPP with 
K users uniformly distributed in the corresponding Poisson-Voronoi cells. 


Both approaches to improve EE are associated with non¬ 
trivial tradeoffs. In this paper, we optimize the deployment of a 
cellular network for higher EE on the uplink. The optimization 
variables are the BS density, the number of antennas and active 
users per BS, the transmit powers, and the pilot reuse. We 
prove the fundamental interplay between these parameters and 
show how to achieve high EE. Some previous works are the 
analytic single-cell analysis in ||7|, 0 and the numerical multi¬ 
cell analysis in 0, 0. This is a continuation of our prior work 
in pO) , which was limited to perfect channel state information. 

II. System Model 

We consider a cellular network that is designed to serve 
users from a heterogeneous user distribution. To this end, we 
adopt the stochastic geometry model in 0, where the BSs 
are distributed in according to a homogeneous PPP d);, of 
intensity A [BSs per km^]. More precisely, this means that in 
any area of size A km^, the number of BSs is a random number 
from a Poisson distribution with mean value XA. These BSs 
are uniformly and independently distributed over the area. 

Each BS has M antennas and serves K single-antenna user 
equipments (UEs). Each UE connects to its closest BS, hence 
the coverage area of a BS is its Poisson-Voronoi cell; see 
Pig. [T] The UEs are assumed to be uniformly distributed in the 
cells. Since there are K UEs per cell, small cells have higher 
user density (per km^) than larger cells, giving a judicious 
deployment based on a heterogenous user distribution!^ 

We consider uplink transmission where the time-frequency 
resources are divided into blocks of Tc seconds and Wc Hz. 
This leaves room for a total number of S' = TcWc transmission 
symbols for pilots and data. The channel response between 
BS/ G <I>A and UE k in cell j is modeled in each block as 

’Alternatively, the UEs could have been distributed according to an inde¬ 
pendent homogeneous PPR This is less sensible since smaller cells would on 
average have fewer UEs than larger cells, which contradicts the main principle 
of densifying networks mainly at places with high user loads. 
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a random Rayleigh fading vector hijk ~ CA/’(0, 

The variance is given by the Euclidean distance dijk [km], the 
pathloss exponent a > 2, and the parameter to that determines 
the propagation loss at the reference distance of 1 km. 

The translation invariance of PPPs allows to perform statis¬ 
tical performance analysis for a typical UE at the origin, which 
is statistically representative for any other UE in the network 
0 - Assume that this typical UE has the arbitrary index k 
and is connected to a typical BS, denoted as BSq S ‘I’a- Then, 
the following basic properties hold (see for example p^): 

Lemma 1. The distance from the typical UE to its serving BS 
is doQk ^ Rayleigh (The BSs of the other cells form a 
homogeneous PPP T'a = $a\{BSo} in {x S : ||x|| >fioofc}- 

The received signal yo G at BSq is modeled as 

K K 

yo = X! V^SOihooi + X! X! + 'lo (1) 

i=l i=l 

where no ~ CN{Q, ct^Im) accounts for the receiver noise. The 
arbitrary symbol transmitted by UE i in cell j is denoted by 
Sji G C and is normalized to unit power (i.e., E{|sjip} = 1), 
while the corresponding transmit power pji > 0 is selected 
as Pji = ptcdfji with p > 0 being a power-control coefficient 
that will be optimized later on. This is a power-control policy 
based on statistical channel-inversion that gives an average 
effective channel gain of E = Mp for all UEs 

(irrespective of their locations), which is necessary to avoid 
near-far blockage in uplink multi-user MIMO systems. 

Pilot transmission is used to estimate channels. We assume 
that there are B pilot symbols per block, where K < B < S. 
Each BS selects K different pilots uniformly at random in 
each block and we call (3 = ^ > 1 the pilot reuse factor. 
Hence, the probability that BSj has a UE (with index k for 
convenience) that reuses the pilot of the typical UE is l/f3. Let 
Xokj G {0,1} be the corresponding random variable, where 
Xokj = 1 has probability 1//3 and means that cell j reuses the 
pilot of the typical UE, then BSq receives the signal 


III. Problem Eormulation 

In this section, we formulate the EE maximization prob¬ 
lem considered in this paper. Since the ergodic capacity of 
a network with imperfect channel knowledge and inter-cell 
interference modeled as a shot-noise process is not known yet 
the first step towards formulating the EE is thus to obtain 
a tractable SE expression. We derive the following lower bound 
using the approach from 0 and Jensen’s inequality. 

Proposition 1. If MRC is employed, a lower bound on the UL 
average SE [bit/symbol/user] is 

(l-^)log 2 (l + SINR) (5) 

where SINR is given by at the top of the page. 

This lower SE bound holds for any /3 > 1 such that f3K < 
S, since the pilot signals need to be contained in a coherence 
block. Observe that f3K does not need to be an integer since 
an arbitrary f3K can be achieved by switching (for appropriate 
fractions of time) between the closest smaller integer and the 
closest larger integer. 

This work aims at maximizing the EE with respect to the 
optimization variables 9 = {(3, p, X, K, M). The area spectral 
efficiency (ASE) [bit/symboEkm^] is computed as 


ASE( 6 >) = AA:(i - ^jlogJl -f SINR ). (7) 

The area energy consumption (AEC) [Joule/symboEkm^] ac¬ 
counts for radiated signal energy, dissipation in analog hard¬ 
ware, digital signal processing, backhaul signaling, and various 
overhead costs (e.g., cooling and power supply losses). A 
detailed model of these factors was given in Q: 
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where we have used that 
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during pilot transmission from the typical UE. The minimum 
mean-squared error (MMSE) estimate of hoofc is 113| 
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and the estimation error Ahoofc = hoofc ~ hoofc is distributed 
as Ahoofc CJV{0, Coofc) with 
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We assume that the BSs use the estimated channels for 
maximum ratio combining (MRC). In particular, the data 
symbols transmitted by the typical UE are detected at BSo by 
correlating the received signal in Q with the MMSE estimate 
as rofc = fc'oofchoofcyo, where i/Qofc G C is a scaling factor 1131. 


^{Pji} = Pw 
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and that each UE transmits S — j3K data symbols and one pilot 
symbol per block p3] . The amplifier efficiency is p G (0,1], 
Co models the static energy consumption of a BS, and 
models the energy consumption of the BS transceiver chains. 
Moreover, CiK+ViMK accounts for the energy consumed by 
a UE and the BS signal processing. The forthcoming analysis 
holds for any positive value of these parameters, but examples 
are given in Table |I] 

The objective of this paper is to solve the following 
constrained EE maximization problem: 


maximize EE( 6 ») = 

see AEC(0) 

subject to SINR > 7 


( 10 ) 


where 0 = {0 : p > 0, A > 0,/? > 1, (M, AT) G Z+, < S} 

is the set of feasible optimization variables and the parameter 






















7 > 0 imposes an average SE constraint of log 2 (l + 7 ) 
[bit/symbol/user], where the average is computed with respect 
to both BS and UE locations. The need for such a constraint 
comes from the observation that an unconstrained EE maxi¬ 
mization may lead to very low SEs per UE. 


IV. Energy-Efficiency Maximization 
In this section, we solve the EE maximization problem 
in The proofs of all theorems are omitted due to space 
limitations but can be found in |[T3]. 


A. Optimal Pilot Reuse Factor j3 

We begin by deriving the optimal value of the pilot reuse 
factor P when the other optimization variables are fixed. As 
stated in the following theorem, the EE-optimal value of /3 is 
such that the SINR constraint is met with equality. 


Theorem 1. Consider any set {p, X, M, K} for which the 
problem is feasible for some (3. The SINR is an increasing 
function of /? and the SINR constraint in is satisfied by 
selecting 


(3* 


Bij 

M — B2y 
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This theorem explains how the EE-optimal pilot reuse 
factor j3* depends on the other system parameters. Recall 
that larger leads to higher estimation accuracy. Erom o- 
([T3J, it is easily seen that, to guarantee a certain SINR, /3* 
increases with K (since the interference increases with more 
UEs) whereas (3* decreases with p as well as with M (as 
it follows from taking the derivative of /3* with respect to 
M). This is because both p and M amplify the desired signal, 
which as consequence improves channel estimation and makes 
the system operate in a less noise limited regime. 


B. Optimal BS Density and Radiated Power 

Based on (3* in Theorem [T] the optimal values for the BS 
density A and the power control coefficient p are as follows. 

Theorem 2. Define p = Xp for p > 0 and consider any set 
of {p, M, K} for which the problem is feasible using j3*. 
The EE is then a monotonically increasing function of X and 
is maximized as A —I 00 . 

This theorem proves that from an EE perspective it is 
preferable to have as high BS density as possible. This might 
sound counterintuitive at first, since the inter-cell interference 
typically grows as the cells become smaller. However, for any 
a > 2 the adopted power control policy reduces the average 
transmit power as with the BS density A (as it 

follows from (|^). Thus, the desired and interfering signals 
are amplified equally much on average; in other words, it is 
only the impact of noise that vanishes as A increases. 

It is clearly unreasonable to have an infinitely high BS 
density, but the numerical results in Section |V] show that the 
asymptotic limit is approached already at A = 10 BS/km^, 
which is a modest number. In ultra-dense deployments, A is 


expected to be one or several orders of magnitude larger than 
this. Hence, we can proceed the EE maximization analysis by 
actually letting A —I 00 without loss of realism. 


C. Optimal Number of Antennas and UEs per BS 

By using Theorems and the original EE maximization 
problem in ( fT0| ) has been reduced tc[^ 


maximize 

M,KGZ+ 
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where we have defined (as obtained from ( [T^ and ( [T3] l when 
p —I 00 ) 
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To find the optimal values for M and K, an integer-relaxed 
version of ( [T4| ) is first considered where M and K can be any 
positive real-valued scalars. The integer-valued solutions are 
extracted from the relaxed problem. Eor analytic tractability, 
we first fix the number of BS antennas per UE, c = M/K, 
and find the EE-maximizing value of K as follows. 

Theorem 3. Consider the optimization problem ([T4J where M 
and K are relaxed to be real-valued variables. For any fixed 
c = MIK > 0 such that the relaxed problem is feasible for 
some K, the EE is maximized by 


K* = 
where 


\J{GCo) + CqDic -f CqG (Cl -f Dqc) — GCq 
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Similarly, if K is fixed, then the EE-maximizing value of 
M is obtained as follows. 


Theorem 4. Consider the optimization problem ( |14| l where M 
and K are relaxed to be real-valued variables. For any fixed 
K > 0 such that the relaxed problem is feasible, the EE is 
maximized by M* in at the top of the next page, if it 
satisfies the constraint in Otherwise, it is maximized by 
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The following parameters are used in these expressions: 
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^ We have also dropped the constraint /3* < S/K since the EE becomes 
negative when this is not satisfied, thus the maximal EE is not affected. 
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Theorems |3] and |4] show how K and M are related at the 
EE-optimal points. We notice that that M* increases with K, 
and vice versa. An important question is whether the EE is 
maximized at AT = 1 (single-user transmission) or at AT > 1 
(multi-user transmission). The answer is ultimately determined 
by the hardware characteristics: Erom ( fTS) ! it is found that 
K* increases with the static energy consumption Co, while 
it decreases with Ci, Vq, and Vi that are the terms of the 
AEC that increase with AT and M. Similarly, M* increases 
with Co and Ci, but decreases with T>o and Vi. The intuition 
is that more hardware should be turned on (i.e., BS antennas 
and UEs per cell) only if the increase in circuit power has a 
marginal effect on the total AEC. Similarly, with a larger static 
consumption we can afford to turn on more BS antennas and 
UEs since the relative power cost is lower. In addition, we note 
that M* is an increasing function of 7 , since deploying more 
antennas is a reliable way to achieve higher rates. 

We can now devise an alternating optimization algorithm 
that solves the integer-relaxed EE maximization problem: 

1) Eind a feasible starting point (M, AT) to (d; _ 

2) Optimize K for a fixed c = M/K using Theorem 

3) Optimize M for a fixed K using Theorem]^ 

4) Repeat 2)-3) until convergence is achieved. 

This algorithm can be shown to converge to the global 
optimal solution {K**,M**) G of the integer-relaxed 
version of the EE maximization problem ( [l4| ); see | [T3) for 
details. It is hard to find an equally strong result for the integer¬ 
valued global optimum to d, but (K**,M**) is a good 
starting point for finding the best integers AT and M. In fact, 
one can show that the EE is quasi-concave with respect to AT 
and M, which implies that the integer-solution is contained in 
a convex level set around (AT**, M**), which effectively limits 
the distance from the solution to the relaxed problem. 

To summarize, the EE maximization problem in ( [T0| l was 
solved by selecting (3 such the SINR constraint is satisfied with 
equality, letting A —> 00 (which removes the impact of p), and 
devising an alternating optimization algorithm that provides 
the real-valued K and M that maximize the EE. The final 
solution is then found by searching through the 

integer points in the vicinity of the real-valued solution, and 
capitalizing on the quasi-concavity. In the process of solving 
Theorems [TJQ have also exposed the fundamental inter¬ 
play between the optimization variables and how they depend 
on the hardware characteristics and propagation parameters. 

V. Numerical Examples 


TABLE I. Simulation Parameters 


Parameter 

Symbol 

Value 

Pathloss exponent 

Ol 

3.76 

Coherence block length 

s 

400 

Propagation loss at 1 km 

LO 

130 dB 

Power amplifier efficiency 

V 

0.39 

Symbol time 

r 

jTt [s/symbol] 

Static energy consumption 

Co 

low ■ r [J/symbol] 

Circuit energy per active user 

Cl 

O.lW-r [J/symbol] 

Circuit energy per BS antenna 

•Do 

IW - r [J/symbol] 

Signal processing coefficient 

Di 

1.56 ■ 10~^*^ [J/symbol] 

Noise variance 


10"^“ [J/symboI] 



Fig. 2. Optimized Energy efficiency (in Mbit/Joule) as a function of the BS 
density for different SINR constraints. 
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In this section, we provide numerical results that illustrate 
and validate the analytic results of previous sections. The 
optimization problem in ( [TOl l depends on a number of parame¬ 
ters that describe the propagation and hardware characteristics. 
Table [^provides the parameter values used in the simulations. 

Fig. 1^ shows the EE as a function of A, when the other 
variables are optimized numerically for each given value of A. 
We show results for the lower bound on the SE in Proposition 
and an upper bound obtained by numerical averaging over 
the BS and UE positions for a finite number of interferers. We 
consider the SINR constraints 7 G {1, 3, 7} and in all three 


Fig. 3. Energy efficiency (in Mbit/Joule) for 7 = 3 . The optimum is star- 
marked and the convergence of the alternating algorithm is shown with circles. 

cases there is only a small gap between the lower and upper 
bounds. The important message from Eig. is that the EE 
increases with A, but begins to saturate around A = 10. Hence, 
EE maximization based on letting A —?► oo gives realistic 
results in both contemporary and future urban deployments. 

Eig. shows the EE for 7 = 3 as a function of M and AT, 
when A, p, and (3 are optimized. The maximal EE is 10.156 
Mbit/Joule and is achieved by (M*, AT*, P*) = (89,10, 7.24). 





























Fig. 4. Energy efficiency (in Mbit/Joule) as a function of the UE density (i. 
The EE is optimized according to with the extra constraint fi = KX, ot 
only with respect to (A, fi, p) for given M and K in the reference cases. 



Fig. 5. Optimized BS density (in BS/km^) as a function of the UE density 
(in UE/km^). The system is optimized in the same way as in Fig. 


This solution is a massive MIMO configuration with a high 
pilot reuse factor. The rationale is that large M and give 
good interference suppression, while the large iT is a way to 
share the energy costs associated with a BS between UEs. 

The alternating optimization algorithm is also illustrated 
in Fig. It is initiated at (M, K) = (20,1) and converges 
in three iterations to the real-valued solution {M**,K**) = 
(91.2,10.2) with an EE of 10.375 Mbit/Joule, which is 0.013% 
higher than the EE at the integer-valued solution. Hence, the 
EE performance is quite flat around the global optimum; it 
decreases quickly with K but slowly with M. 

A. EE Maximization for a Given UE Density 

Next, we study the tradeoff between massive MIMO and 
small cells when a cellular network is deployed to cover a 
heterogeneous UE density of p UE/km^. Mathematically, this 
amounts to solving ( [T0| ) with the additional constraint fi = K\, 
which can be easily solved numerically. We consider the range 
e [10^,10®] since the METIS project report predicts 
future UE densities from p = 10^ UE/km^ (in rural areas) to 
/i = 10^ UE/km^ (in shopping malls). 

Fig. 0 shows the EE as a function of the UE density /r 
for the average SINR level 7 = 3, while Fig. shows the 
corresponding BS density. Apart from the optimal solution, we 
consider two reference cases: single-user single-input multiple- 
output (SIMO) transmission with (M^K) = (10,1); and 
massive MIMO with {M,K) = (89,10). 


We notice that the EE level is roughly independent of the 
UE density for p, > 100, which includes all the cases predicted 
in 0. Single-user transmission is attractive for very low UE 
densities, while massive MIMO transmission is the best choice 
in the practical cases: Fig. shows how it improves the EE 
by a factor 3 and Fig. shows that for a given UE density 
we can have 10 times fewer BSs per km^. The latter has the 
benefit of greatly reducing the deployment costs. 

VI. Conclusion 

Massive MIMO and small cells are two approaches to 
improve the energy efficiency of wireless networks. The main 
benefit of small cells is the reduced propagation losses, while 
the main benefit of massive MIMO is the interference suppres¬ 
sion among the UEs that share the energy costs associated with 
the serving BS. Hence, these techniques are not competing 
but have complementary features that should be combined to 
achieve the maximal EE in future wireless networks; massive 
MIMO and small cells are a perfect match for marriage. 
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